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Abstract 

We study several anisotropic inflationary models and their implications 
for the observed violation of statistical isotropy in the CMBR data. In two 
of these models the anisotropy decays very quickly during the inflationary 
phase of expansion. We explicitly show that these models lead to violation 
of isotropy only for low / CMBR modes. Our primary aim is to flt the 
observed alignment of / = 2, 3 multipoles to the theoretical models. We 
use two measures, based on the power tensor, which contains information 
about the alignment of each multipole, to quantify the anisotropy in data. 
One of the measures uses the dispersion in eigenvalues of the power tensor. 
We also deflne another measure which tests the overall correlation between 
two different multipoles. We perturbatively compute these measures of 
anisotropy and fix the theoretical parameters by making a best flt to 
Z = 2, 3 multipoles. We show that some of the models studied are able to 
consistently explain the observed violation of statistical isotropy. 



1 Introduction 

There exist many observations which suggest violation of the cosmological prin- 
ciple. These include dipole anisotropy in radio polarizations [lj|2], large scale 
alignment of optical polarizations [S], alignment oi I = 2,3 multipoles of the 
cosmic microwave background radiation |3HS] i dipole anisotropy in distribution 
of radio galaxies [7j and cluster peculiar velocities [S] . Remarkably all of these 
indicate a preferred direction towards the Virgo cluster, close to the observed 
CMB dipole. This coincidence of diverse axes was first pointed out in [9l[T0]. 

There also exit several other indications of violations of isotropy. These in- 
clude dipole modulations to CMB including north-south ecliptic power asymme- 
try prUTT) , quadrupolar power modulation to CMB |16ll81l^ , dipole anisotropy 
in galaxy distribution [5T], and parity violation in the CMB data and 
in the handedness of the spiral galaxies [24]. The indications for a preferred 
direction in diverse data sets have motivated many theoretical studies of infla- 
tionary models which violate statistical isotropy and homogeneity, giving rise to 
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a direction dependent power spectrum [T51[^5H5S] . Anisotropic inflation is likely 
to lead to observable effects in the CMBR j37j . 

It is well accepted that the observed Universe is neither isotropic nor ho- 
mogeneous at early times. It may evolve into a de Sitter space as it expands 
in presence of a positive cosmological constant. This has been explicitly es- 
tablished for Bianchi models, which represent a homogeneous and anisotropic 
Universe. An anisotropic Bianchi model is expected to become isotropic within 

a time scale of order during the early stages of inflation [35] . It was pointed 
out in [33], that even this brief early phase of anisotropic expansion can lead 
to violations of statistical isotropy in the current era. It was shown that for a 
wide range of parameters, modes generated during the anisotropic inflationary 
period can re-enter the horizon at recent times and hence affect observations. 
For some parameter range they may re-enter even before decoupling of radiation 
from matter [39] . In the current paper we shall show that the anisotropic modes 
arising during the very early phase of inflation can consistently explain the align- 
ment of CMB quadrupole and octopole. We point out that our explanation is 
within the framework of the inflationary Big Bang cosmology. 

There exist several techniques to study the violation of statistical isotropy 
in the CMB data. These include bipolar spherical harmonics [iUlBTj , the power 
tensor [HISS], multipole vectors :5j etc. Here we shall use the power tensor, 
which is defined by, 

^'^■^^^ " ITTTT) ^ (/™|J,|Zm')(/m"|JjlMa/m'aL" (1) 

m,m',m" 

where aim are the standard harmonic coefficients of the temperature anisotropy, 
AT(ji)/T, and Ji are the angular momentum operators in spin— ^ representation. 
This tensor is a second rank real symmetric matrix. Its eigenvectors define 
an invariant frame for each multipole and the corresponding eigenvalues yield 
the power along each of these vectors. For a given multipole, the eigenvector 
associated with the largest eigenvalue of the power tensor is called the principal 
eigenvector (PEV). This defines the preferred direction for that multipole [HI 
HI] . Statistical isotropy implies that the two point correlations of the harmonic 
coefficients obey, 

(a/„ia;*„j,) = 5,;'(5„„/C/ (2) 

where C/ is the power spectrum and the angular brackets denote an ensemble 
average. For the power tensor, this implies, 

(^..(0) = y%- (3) 
In the case of an anisotropic background metric, 

{aimOiiyn') = Cu'mm' , (4) 

i.e. the correlations are not diagonal as in the isotropic case. This will im- 
ply that, in contrast to Eq. (3] the power tensor is not proportional to Sij. 
Furthermore it will also generate correlations in the power tensor for different 
multipoles. 

The alignment of the CMB quadruple and octopole can be tested in the 
power tensor approach by determining the corresponding PEVs ^1421143] . These 
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are found to be aligned very closely with one another. In the 7 year WMAP 
data release, the significance of alignment is found to be better than 4tT CL in 
the ILC map [13]. The multipoles, 1 — 2,3, by themselves, however, do not give 
any indication of violation of statistical isotropy. The hot and cold spots, both 
for / = 2 and / = 3 do lie roughly in a single plane. However the effect is not 
statistically significant. It is only the alignment of the two multipoles which is 
significant. 

In our theoretical approach we work in a preferred coordinate system with 
the preferred axis aligned along the PEV corresponding to / = 2. We use two 
different approaches to test for alignment. In one approach we determine the 
PEV for each multipole and then test for its alignment. In the second approach 
we follow a blind procedure which directly tests the overlap of the two power 
tensors. The overlap is defined as the trace of the matrix {A'^{l)A{l'))^j. We 
may express this ensemble average as. 



^ I I' 

" N(l]N(l') ^ ^ ^ {lmi\Jk\lm2)*{lm3\J,\lmi)* 
{l'm\\Jk\l'm'2){l'm'3\Jj\l'm{){a*i^^ai^.,ai,^>yi,^,^) . (5) 

where the normalization N{1) = l{l + 1). Here on the right hand side the 
subscript n takes values 1,2,3, i.e. the sums are over mi, m2, ma, m'j^, 7113. 
The corresponding statistic is defined as, 

Su.^^Jr{A+{l)A{l')), (6) 

For the isotropic case, the expectation value on the right hand side factorizes 
and we obtain, 

S*;;/ (isotropic) = 1 . (7) 
Here we shall only be interested in the case V = 1 + 1 and hence use the statistic, 

Si = j^Tr{A+{l)A{l + l)). (8) 

The blind approach has the advantage that it tests the overall correlation of 
the two multipoles rather than testing a specific feature. Hence it does not 
suffer from the criticism of being an a posteriori statistic [2^. As expected, this 
statistic leads to a smaller significance of alignment between ^ = 2, 3. 

The power spectrum P{k) for the primordial density perturbation 5{k) is 
defined as 

{5{k)5*{k')) =P{k)5'^{k~k') (9) 
where 5{k) is the Fourier amplitude of the density inhomogeneities, 

^^-WY. h^^y~'"'d'k. (10) 
p[x) [2tiY J 

In the case of a homogeneous and isotropic FRW metric, the fluctuations are 
statistically isotropic and the primordial power spectrum P{k) depends only on 
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the magnitude of the wave vector k. However, in the case of anisotropic models, 
which we shall consider in this paper, it will also depend on the direction. The 
direction dependent contribution to -P(fc) has been obtained in a particular 
anisotropic model in Ref. [15] . Here we examine this model as well as two other 
models of anisotropic inflation. 



2 Anisotropic Metric 

We consider three different anisotropic metrics. The axis of anisotropy is taken 
to be the z— axis. The flrst model may be expressed as, 

ds^ ^ dt^ -2y/^dzdt- a^{t){dx^ + dy^ + dz^) Model I (11) 

where a{t) is the scale factor and a is the anisotropic parameter which controls 
deviation from statistical isotropy. The second model we consider is same as 
the one used in jl8| . The corresponding metric is given by, 

ds^ = df - a^{t)dxl-b^{t)dz^ Model H (12) 

where dx'\^ = dx^ + dy^. This model is anisotropic throughout inflation. The 
third model is a modiflcation of this model such that the anisotropic term be- 
comes insignificant after one cfold. The line elements may be expressed as, 

ds^ ^ dt^ - al{t){dx^ + dy^) - al{t)dz^ Modellll (13) 

where 02 (t) = ai{t) + a, a being a constant, independent of time. We also 
define, d{t) — {al(t)a2(t))'s . 

Here we shall present the calculation of the alignment between multipoles 
for the case of model I in detail. The calculation for model II and HI proceeds 
along the same line and for these models we only give the final results. We shall 
treat the anisotropic term perturbatively, following the treatment of [IH]. For 
models I and HI this term is expected to be significant only at very early times. 
For model II, this term would contribute during the entire inflationary period. 
In all the models, the anisotropy parameter may be fixed by fitting the data for 
1 — 2,5 multipoles. The reliability of perturbation theory will depend on the 
magnitude of this parameter. 

The components of the Einstein tensor for the perturbed metric correspond- 
ing to model I are 

_ a{t)[5<T + 3a^{t)]d^{t)-2<T[a + a^{t)d{t)] 

Goo — 7777 — ; — 277yv2 

a{t){a + a'^{t)Y 



^ [3a + a'{t)]a\t) + 2a{t)[a + a\t)\d{t) 

<-^03 — -yo^ F ; TTTTT^ 



2 [3a + a'{t)]d\t) + 2a{t)[<j + a\t)\d{t) 
[fj -I- a"^(t)J"^ 

G30 — Go3 (17) 
Gil — G22 = G33 . (18) 

The perfect fluid energy momentum tensor is deflned as, 

T^^u = (/O + P)Uf^u^ - . (19) 
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This leads to the fohowing off-diagonal terms, 

To3=T3o = -pV^- (20) 

The Einstein equations, to zeroth order in the anisotropic parameter y/a, are 
given by, 

-V + 2- = -SttGP (21) 
a J a 

3 Q = SttGp (22) 

where the dot refers to the derivative with respect to the cosmic time t. Hence 
with P = —p we get standard inflationary solution, 

a{t)^aie"* (23) 

where H is the Hubble constant during inflation. 

3 Quantization of field and power spectrum 

In this section we discuss the quantization of the massless scalar field '44' and 
obtain the power spectrum perturbativcly to first order in the anisotropy pa- 
rameter, y^. The action may be written as, 

S^,^ J d'^x C (24) 

where the Lagrangian density, 

£ = lV^9'"'d,cj,d,<j>. (25) 

The scalar field (t>{x, t) is quantized in the standard manner by expanding in 
Fourier modes, 

Mx,t) = / (0 {e''-'Mt)ak + e-^'-UUM) ■ (26) 

The creation and annihilation operators and satisfy the commutation 
relation [a^, a^,] — {2TT)^d{k — k'). The Euler-Lagrange equation of motion, 

^ a.(g''''y^a^^) = 0, (27) 



V-5 

for the isotropic FRW metric reduces to 



2 

3i70-^0 = O. (28) 



In conformal time, this can be written as. 



v'i+[k'-^^]v,^Q, (29) 
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where, Vk — a<j)k and ' — denotes the derivative with respect to the conformal 
time. The solution to the unperturbed equations, (f>1{rj), is given by, 



H 



rj j exp(— ifcry) . 



The Hamiltonian in interaction-picture is found to be 

' d(j)i{x,vi) d(j)i{x,ri) 



Hi = J d^x 



drj dz 

The two point correlations, to first order, are given by, [T51H5] . 

{(j){xi,t)(l)[x2,t)) = {(I)i{xi,t)(j)i{x2,t)) + 

i f dt' {[Hi{t'),(i,i{xi,t)d^i{x2,t)]) 



(30) 



(31) 



(32) 



This leads to, 

{(j){xi,t)(f){x2,t)) 

where 
and 

AP(fc) = -2V^ 
Here we have used, 



(27r)= 



■ exp 



ik-[xi-x2) [P.so(fc) + fc,AP(fc)] (33) 



PUk)c,\4>f{r^)\ 



a{ri')dri' 



V = 



dt 



-Ht 



a{t) Haj 
The complete expression for AP{k) is given by, 



AP(fc) = P,,„(fc) ( ^ ) <! Hrji2 + k^rj^' 



(34) 
(35) 

(36) 



ai 



-1 - iajHri + k^rf) cos 



2k 



aiH 



2a/fc 



[Ak'^r] + 3a/i7(-l + fc^ry^)] sin 



2k 



aiH 



.(37) 



Here we assume that the universe was anisotropic at early times. The perturba- 
tions which are generated during the early anisotropic phase may re-enter the 
horizon during matter or radiation dominated phases. These modes crossed the 
horizon very early during inflation. This implies k ^ ajH. Taking the limit 
|fc?7| << f, we find that, 



Piso {k) 



2fc3 



2k 



f 2k \ f3Hai\ . 



(38) 
(39) 
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Finally, the modified direction dependent power spectrum is given by, 

P'{k) = P,soik)[l + ik- z)g{k)] Model I (40) 

where k is the unit vector along the direction of fc, z is the unit vector along 
the preferred direction and the function g{k) is defined as, 



9ik) 



ai 



\aiH J \ 2k J \aiH J 



(41) 



The corresponding results in model II have already been obtained in Ref. 
The power spectrum can be expressed as, 



P\k) = P,so{k) l+g{k)Ck-nY 



Model II 



(42) 



where fi — (0,0, 1) in the notation of Ref. [18j. The function, g{k), at the end 
of inflation (time i*) is given by, 
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9ik) = -e/flog 



qjt*) 

H 



(43) 



H = ^{2Ha + Hb) and the physical wavelength q{t^) ~ jrri- "^^^ anisotropic 



parameter is defined as 



2 I Hi, — Ha 



(44) 



Here we have neglected the perturbative correction to the isotropic power spec- 
trum. 

In the case of model III, the Hamiltonian in interaction picture is 



Hj{t) 



1 



d-'x- { (a(t) - a2{t)) 



a{t) 



dx 



dy 



al{t)\ f d(f> 



dz 



(45) 



The two point correlation function is found to be 

d^k 



{(I){xi,t)(j){x2,t)) = 

where 



(27r) 



■ exp 



^k■{x^-X2) [PUk) + kfAPik)] (46) 

(47) 



AP(fc)^-2a/ ^ a{r^')drj'[2Im(^cf>lirj'f4*{r,) 

The modified power spectrum is given by, 

P'{k)^Pi^{k) [l + {k-zfg{k) 

where, 



(48) 



with 

Si(x)^ fdx'^. (50) 
Jo ^ 

Here the isotropic power spectrum, P!gg{k), includes a perturbative correction 
to Piso{k). It is given by, 



(51) 



4 Effect of the power spectrum on the Cosmic 
Microwave Background 

The temperature perturbation can be decomposed in terms of spherical har- 
monics, Yimiv), as follows, 

5T °° ' 

— (f,p,?7) = ^ ^ a;„yi„(p) . (52) 

l—l m— — l 

This temperature anisotropy ^ can be related to the primordial density fluc- 
tuations as : 

5T 



T 



{p) = J dkY,^^{-iypi{k-mmi{k) (53) 



where Pi is the Legendre polynomial of order I and Qi (fc) is the transfer function 
which correlates initial fluctuations to the observed temperature anisotropies. 
The expansion coefficients aim can be calculated using. 



air. 



J dnY;^ip)^-fix,p,r,). (54) 



Wc now compute the two point correlation function of aim's using the directional 
dependent power spectrum. We have, 

{aima'l^m'} = {aimat,m')iso + {(^lma*l'm')aniso (55) 

where the first term corresponds to the isotropic case. It is given by 

/•oo 

{aimalm')iso = ^ll'^mm' I k'^ dkPiso{k)&f (k) . (56) 

Jo 

The second term which contains departures from statistical isotropy is given by, 

{aima:,m')aniso = i-it''ilm;l'^' / k^ dkPiso{k)g{k)ei{k)ei, {k) (57) 

Jo 

where we have used the spherical components of the unit vector n as 

n+ = -(— -^), n_ = ( — no = • (58) 
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The constant iimd'm'i which encodes the correlation between different multi- 
poles, is defined as 



/ 47r 



where, 



\l-m + l){l-m + 2) 
(2/ + l)(2/ + 3) 



' {l+m-l){l + m) 
(2Z + 1)(2/-1) 



5v. 



i-i 



'{1 + in + + t,i) 
(2Z + l)(2Z + 3) 



Si' 



'{l-m-l)il-m) 

(2/ + l)(2/-l) 



Si'.i-i 



Sm' . 



' {I - m+l){l + m+l) 
(2Z + l)(2Z + 3) 



(5r 



' (I — m){l + m) 
(2Z + 1)(2Z- 1) 



Similarly, using = (— we obtain. 



Jo 



and 



f 

Jo 



(59) 



(60) 



(61) 



(62) 



k^dkPiso{k)ef{k) (63) 



edkPiso{k)g{k)Qi{k)Qi,{k) . 



We express ^im;i'-m' as, 



lm:l' —m' ' 



(64) 
(65) 
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where, 



^m'. — (m— 1) 



(/-m + l)(?-m + 2) 
(2Z + l)(2Z + 3) 



'(/ + m- l)(Z + m) 
(2? + 1)(2/- 1) 



(66) 



2 ^m' , — (m+l) 



'(? + m + l)(/ + m + 2) 

(2; + i)(2; + 3) 



'{I - m - - m) 
{2l + l){2l- 1) 



(67) 



'{l-m+ + m + 1) 
(2^ + l)(2Z + 3) 



[I -m){l + m) 



Si', 



/-I 



y (2; + i)(2;- 1) 

Computing the trace of {A'^{l)A{l + 1)), we find 

1 



(68) 



Tr{A+il)Ail + l)) 



4Z(/ + l)2(/ + 2) 



E 



d+(l,mi)d {I + l,m5)Sm2.rni + lS„u 



ms — 1 



nil ■:f^i-2 ,m5 



d'^{l,m3)d {l + l,m4)SrnuTn3 + lSme.Tni-l+d (^ TO3)d+ (Z + 1, 7714) 



^ml,m3 — l'^m6, 7714 + 1 



^7712 ,7713 ^7715 , 7716 



edkP,,o{k)ef{k) / k^dkP,so{k)&f+i{k) 



+ (Cim2;(i+l)meCm3;(i+l)m5 + ?im2;(i+l)"iBCm3;(i+l)m6) 



edkP,so{k)g{k)Qi{k)Qi+^{k) 



(69) 



where we defined (i^(/, m) = ^(/ — m){l + m + 1) and (i^(/, m) = {I + m){l — m + 1). 
Note that the sum in the above equation is over nii, i — 1, 2, 6. Here we take 
the transfer function for large scales as &i{k) — ^Ji[k{rio — rjd)], where 770 is 
the conformal time today and rjd is the conformal time at last scattering. Since 
Vd ^ Vo we take Ji[k{riQ ~ 77^)] « J/[fc?7o]- It is convenient to define, x — krjQ. In 
terms of x the function g(k) in Eq. 1411 mav be expressed as, 



g{x) = 

ai 



f 2x \ iH-qaai . f 2x \ 
\rioajH J 2x \rioaiH J 



(70) 
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The integral in Eq. 15^ involves an oscillatory functions of x. It is clear that 
for large I the integral will get contributions from relatively large values of k. 
For sufficiently large k the integral will be damped due to oscillations in g{x). 
Hence the correlations between different I are expected to be small for large I. 

4.1 Models II and III 

For models II and III the constant S,im-i'm' is the same as in [18 and is given by 
£.lm-l'm' = ^ / dVLk{Yi.^{k))*Yi,m'{k) 



3 



2 



X {n+Yl{k)+n^Y^^{k)+noY^{k)y 

+ 2n_noC;i'™' + "oCn;Z'm' • (71) 

For I = V the contributing terms of these coefficients are 



_ _ ^{P-{m + lY){l + m + 2){l-m) 

Wm' - ^m',m+2 (2Z + 3)(2Z-1) ^ ^ 

ilmil'm' — ivm'dm (73) 
lihn-l'm' - <^rn',m ^21 - 1){21 + 3) ^ ' 

_ {2m+ 1)^(1 + m+l)il-m) 

l^l^-Mm' - ^'^rn'.rn+l (21 ~ 1){21 + 3) ^ ' 

ihn;l'rn' — '"^I'm'-lm (76) 

,00 _ . (2/^ + 2?-2m^-l) 

Following [m we set, | log | ^ 60. Hence, g{k) = 5* = fe// x 60. The 

ensemble average of the power tensor is given by, 

= —1- J2 {lm\Mlm'){lm"\J,\lm) 



ip'lni' ^Ira"^ iso H" {p'lm' ^^lui"^ aniso 



(78) 



In this model we find that {Aij{l)) 7^ In the preferred coordinate system 

we obtain a diagonal matrix with the component corresponding to the preferred 
direction different from the other two. 



5 Results 

We first give the numerical results for model I. In this case we use the statistic 
Si (Eq. [S]), in order to test the correlation between different multipoles. As 
explained in the introduction, this tests the overall correlation of the two mul- 
tipoles rather than testing a particular aspect, such as the alignment of their 
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principal eigenvectors. For an isotropic model, this statistic takes the value 
unity. We have explicitly verified this by using randomly generated CMB data. 
Furthermore the anisotropic term gives zero contribution to Ci in this model. 
The value of this statistic for the observed CMB data is given in the second 
column in Table 1. Here the results are presented for seven year WMAP data. 
We have checked that the nine year data also gives nearly identical results. The 
fourth column in the table gives the significance or P-value, i.e. the probability 
that the correlation seen in data arises as a random fluctuation in a statistically 
isotropic sample. This is computed by comparing the data statistic with 4000 
random realizations of CMBR based on the ACDM model. We find that this 
statistic also indicates a significant correlation between I = 2 and Z = 3 multi- 
poles but doesn't lead to significant result for other multipoles. This agrees with 
earlier results which show a very strong alignment of the principal eigenvectors 
for Z = 2, 3 [13] and a relatively mild effect for remaining multipoles [12]. 

The theoretically computed values of 5*; for model I are given in the third 
column of Table 1. Here we have taken the anisotropy parameter — ^ — 
1.73 in order to fit the observed 5; for / = 2,3. We also set — ^—rr = k- If we 
choose a smaller value for this parameter, we shall obtain larger correlations 
for larger values of /. This choice of parameters fits the observations, as shown 
in Table 1. Our theoretical results match the data for {1,1 + 1) > (3,4) in the 
sense that neither the observed nor the theoretical values indicate significant 
correlation. We do not expect to obtain an exact match since the observed 
values are statistical in nature. Assuming that statistical isotropy is valid for 
I > 4, we expect that each of these multipoles would be aligned randomly and 
hence their correlations are expected to be small and take random values. 



(1,1+1) Si{Data — value) Theoretical — value P ~ value 



(2,3) 


1.487 


1.487 


0.005 


(3,4) 


0.696 


1.279 


0.969 


(4,5) 


1.030 


1.127 


0.416 


(5,6) 


1.100 


1.039 


0.198 


(6,7) 


1.139 


1.004 


0.102 


(7,8) 


1.055 


1.002 


0.270 


(8,9) 


1.007 


1.013 


0.458 


(9,10) 


1.050 


1.021 


0.227 


(10,11) 


0.989 


1.019 


0.568 


(11,12) 


0.934 


1.001 


0.849 


(12,13) 


0.886 


1.003 


0.968 



Table 1: Comparison of the observed and theoretical results for model I. In 
column 2 and 3 we show the observed and theoretical values of the statistic Si . 
The significance or P-value of the observed statistic is shown in column 4. 

We point out that the contribution due to the anisotropic term is relatively 
small compared to unity, obtained with an isotropic model. Hence we expect 
perturbation theory to be qualitatively reliable. The maximum contribution is 
found to be 0.487 for the case of I = (2, 3). We expect higher order corrections 
of the order of (0.487)^, roughly about 25%. 

We next describe the results corresponding to model II. In this case the 
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metric is anisotropic during the entire inflationary period. Hence we expect 
that ah the modes, independent of the values of k, would violate isotropy and 
lead to alignment among CMB multipoles of arbitrary I values. In this case 
we compute the principal eigenvectors of the power tensor of each multipole. 
For each / we expect two equal eigenvalues. The third eigenvalue is found 
to be larger than these two eigenvalues for negative values of the anisotropy 
parameter en- Hence, in this case, the principal eigenvector would align with 
the preferred (z) axis. We fix the value of en by demanding that the theoretical 
and observed values for dispersion in the eigenvalues match with one another. 
A useful measure for dispersion is the power entropy J42^. Let Xa, a = 1,2,3 
represent the three eigenvalues of the power tensor. As discussed in 42_, these 
eigenvalues are positive. We define the normalized eigenvalues 



Sp{l) can take values in the range to log 3, with being the value for maxi- 
mal dispersion and, hence, anisotropy. In the CMB data sample, the value of 
Sp{l) for any I deviates from the perfectly isotropic limit of log 3 due to cosmic 
variance. The cosmic variance becomes small for large I. Hence, as expected, 
the observed value of Sp{l) approaches log 3 for large I We compute the 

power entropy expected for an isotropic sample by using 4000 randomly gener- 
ated CMB data samples. Let Sp{l,data) and Spil^mean) denote, respectively, 
the power entropy of the observed CMB data and the mean of random sam- 
ples. We attribute the difference ASp{l) = Sp{l,data) — Sp{l,mean) to the 
contribution due to the anisotropic term. The theoretical estimate of ASp{l) 
is, therefore, equal to ASp^th{l) = Sp^th — log3. Here Sp^th is computed using 
the mean power tensor given in Eq. [75] and we subtract log 3 since that is the 
value obtained in the limit of perfect isotropy. 

For Z = 2,3 we find that, ASp{2) = -0.0056 and A5p(3) = -0.064 using 
the 9 year WMAP data. Similar results are obtained for 7 year data. For 
I — 2 the observed value is very close to that expected in an isotropic model. 
For / — 3 the deviation is larger but still less than la. Hence the dispersion 
in eigenvalues for the observed data is not statistically significant for Z = 2 or 
I = 3. We nevertheless fix the value of anisotropy parameters in model II so as 
to obtain a best fit to ASp for Z = 2, 3 by minimizing the square error. In future 
it might be more appropriate to make a global fit to the CMBR data in order 
to fix the parameters of this model. The best fit is obtained for cr — —0.0054, 
which leads to ASp^th — —0.034 and root mean square error equal to 0.041. The 
value of ASp,th is found to be independent of I since the anisotropy is equally 
effective for all the multipoles. The observed data, however, shows dominant 
alignment only for I = 2,3 multipoles. Hence we need some mechanism to 
suppress the anisotropic contribution in the remaining multipoles. This can be 
accomplished if the Universe quickly evolves into an isotropic de Sitter phase 
during inflation [39]. This mechanism is implemented in model HI. 

In model HI, ASp^th oscillates, with amplitude decreasing with /. We at- 
tempt to fit the power entropy for Z = 2, 3 in this model by fixing the parameter 
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The power entropy, Sp{l), may be expressed as. 
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Figure 1: The power entropy as a function of I for model III. The results are 
shown for three different choices of parameters, ct' and q — 2/{rioaiH). 



q = 2/(j]QaiH) and determining cr' = a/aj by minimizing square error. How- 
ever we find that the best fit is driven to relatively large values of cr', where the 
perturbation theory becomes unreliable. The largest perturbative contribution 
is obtained for the isotropic part of the power spectrum. Hence here we only 
give results by fixing the parameter a' = 1.0. The results for Sp^th as a function 
of I for q — 0.8, 1.0, 1.2 with a' = 1.0 arc shown in Fig. [TJ The figure clearly 
shows the oscillations in Sp^th- We also find that it approaches log 3 for large 
Z, as expected. The root mean square error for cr' = 1 and q — 0.8, 1.0, 1.2 is 
found to be 0.058, 0.056 and 0.055 respectively. We postpone a more detailed 
fit, involving larger values of I and perhaps polarization data in order to obtain 
a global minima to future research. In Fig. [T] we notice that the anisotropic 
contribution becomes negligible beyond I ~ 10. This might be related to the 
observations that low multipoles lying in the range 2 < Z < 11, and not just 
I = 2, 3, show alignment with the quadrupole 142 



6 Summary and Conclusions 

We have studied three models of anisotropic inflation in order to explain the 
observed violation of isotropy in CMBR data. In models I and HI, the Universe 
is anisotropic only during the very early stages of inflation and quickly evolves 
into a de Sitter space-time. In model II, the anisotropy is present throughout the 
period of inflationary expansion. In all these models we choose the preferred 
direction same as the principal axis of the CMBR I = 2 mode. This axis 
points roughly in the direction of the Virgo supercluster. Model I leads to 
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direct correlations between multipoles I and ^ + 1 . Hence it leads to alignment 
between ^ = 2, 3. In models II and III, there is no direct correlation between I 
and I + 1 modes. However each multipole is affected such that its principal axis 
is aligned with the preferred axis of the model. Hence these models also lead to 
alignment of the principle axis corresponding to ? = 2, 3 multipoles. In models 
I and HI, the anisotropic contribution affects dominantly the modes with low 
I. For higher / values, this contribution decays rapidly. Model II, in contrast, 
gives equal contribution to all I values. 

We treat the anisotropic term perturbatively in all the three models. We 
compute the power spectrum and measures of alignment of Z — 2,3 multipoles 
at first order in perturbation theory. The measures of alignment use the concept 
of the power tensor [SlIU] . The three eigenvectors of this tensor define a frame 
in real space for each I. The eigenvector corresponding to the largest eigenvalue 
defines the principal axis. In models II and HI the principal axis of both / = 2, 3 
are aligned with the chosen preferred axis in these models and hence are aligned 
with one another. In model II we fix the theoretical parameters by fitting the 
dispersion in eigenvalues in ^ = 2, 3 modes. In model I, the anisotropic term does 
not contribute to the power tensor for any individual value of I . However it leads 
to a correlation between I and l + l modes. We define a new statistic, which tests 
the overall correlation of the power tensor between two adjacent modes. The 
theoretical estimate of this statistic is fitted to the observed data for I = 2,3 
modes. We find that models I and HI provide a better description of data 
since they capture the important feature that the anisotropy decays rapidly for 
large I. We also find that perturbation theory is qualitatively reliable in models 
I and II. The contribution due to higher order terms in model I is likely to 
introduce a correction of about 25%. In model HI, the best fit to the eigenvalue 
dispersion for Z = 2, 3 modes is driven to relatively large values of the anisotropy 
parameter, making perturbation theory unreliable. In any case this might not 
be the best method to fix the model parameters since the dispersion observed 
in data is not statistically significant for these modes. Hence it may be better 
to fix the parameters of both models II and HI by making a global fit to the 
CMBR temperature and polarization data. It may also be useful to perform a 
non-perturbative analysis of the anisotropic metrics. 

The observed data has shown several signals which indicate a preferred di- 
rection, pointing roughly towards the Virgo supercluster. The fact that several 
diverse data sets, including radio, CMBR and optical, lead to the same preferred 
direction gives us some confidence that the observed anisotropy arises due to 
a physical effect. We have related these observations to an early anisotropic 
phase of infiation. It would clearly be of interest to test our models further 
by studying their implications for higher multipoles, CMBR polarization, radio 
and optical polarizations and the large scale structure of Universe. 
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